Accuracy of slow-roll formulae for inflationary 
perturbations: implications for primordial black 
hole formation. 

Sirichai Chongchitnan and George Efstathiou 

Institute of Astronomy. Madingley Road, Cambridge, CBS OHA. United Kingdom. 
E-mail: sc427@cam.ac.uk 

Abstract. 

We investigate the accuracy of the slow-roll approximation for calculating 
perturbation spectra generated during inflation. The Hamilton- Jacobi formalism is 
used to evolve inflationary models with different histories. Models are identified 
for which the scalar power spectra V-ji computed using the Stewart-Lyth slow- 
roll approximation differ from exact numerical calculations using the Mukhanov 
perturbation equation. We then revisit the problem of primordial black holes generated 
by inflation. Hybrid-type inflationary models, in which the inflaton is trapped in the 
minimum of a potential, can produce blue power spectra and an observable abundance 
of primordial black holes. However, this type of model can now be firmly excluded 
from observational constraints on the scalar spectral index on cosmological scales. 
We argue that significant primordial black hole formation in simple inflation models 
requires contrived potentials in which there is a period of fast roll towards the end of 
inflation. For this type of model, the Stewart-Lyth formalism breaks down. Examples 
of such inflationary models and numerical computations of their scalar fluctuation 
spectra are presented. 



PACS numbers: 98.80.Cq 
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1. Introduction 

The simplest models of inflation predict an almost scale-invariant and Gaussian 
distribution of primordial density perturbations as a consequence of the inflaton field 
roUing slowly down a potential of a specific form [1, 2]. Provided that the roUing is 
sufficiently slow, a set of 'slow-roll' approximations allows the primordial perturbations 
to be computed easily. 

However, there is no compelling reason why inflation should be of the slow-roll 
type. Inflationary scenarios that violate the slow-roll conditions have been studied by 
many authors [3, 4, 5, 6, 7]. Furthermore, observations that constrain the amplitude 
of primordial perturbations only probe a relatively small range of scales. It is therefore 
possible that inflation might have been more unusual than that anticipated in the 
simplest theories, and that there may have been a period of fast roll before inflation 
ended. 

In this paper, we analysed inflationary perturbations numerically for single-fleld 
inflationary models with generic evolutionary histories, including cases where slow- 
roll conditions are badly violated. We then compare the numerical results with those 
calculated using slow-roll formulae, and identified cases where the two approaches differ. 

We revisit the problem of primordial black holes generated by inflation. We 
argue that primordial black hole formation in slow-roll inflation models is excluded 
observationally. It is, however, possible to form primordial black holes if there was 
a period of fast roll towards the end of inflation. For this type of model, numerical 
computation of the fluctuations is essential to calculate the abundance of primordial 
black holes. In the flnal section, we give an explicit example of a model which 
satisfles observational constraints on CMB scales, but produces primordial black holes 
in signiflcant abundance. The model is, of course, contrived, but gives an indication of 
the level of flne-tuning needed to produce primordial black holes in single-fleld inflation. 



2. Evolution of scalar perturbations 
2.1. Generating inflationary models 

Given our lack of knowledge of the fundamental physics underlying inflation, we shall 
investigate scalar perturbations in a large number of inflationary models numerically. 
Our approach is based on the Hamilton- Jacobi formalism introduced in [8]. Tensor 
perturbations were treated using this formalism in our earlier work [9] . 

In this approach, the dynamics of inflaton is determined by the Hubble parameter 
which is related to the potential V{(f)) via the Hamilton-Jacobi equation 

iH'mf - ^^H\4>) = -^^(0) . (1) 

The functional form of H{(f)) can be specifled by a hierarchy of 'flow' parameters: 

mh fH'V ml, fH" 
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, (7 = 2?7 - 4e , (2) 

where the definition for the parameter a is motivated by the expression for the scalar 
spectral index [10] 

- 1 ~ (7 - (5 - 3C)e2 - ^(3 - 5C)ae + ^(3 - C)C, (3) 

with C = 4(ln 2 + 7) — 5 ~ 0.0814514. Thus, to lowest order, a describes the departure 
from scale-invariance. 

The flow variables satisfy the following equations (see e.g. [10] and references 
therein) 

1^ =e{a + 2e), 



dN 
da 

dN 

dN 



- e(5(7 + 12e) + 2^ , (4) 
- 1 



a+{£-2)e 'Xh + '^^\h , > 2) 



2 

where N is the number of e-foldings before infiation ends. These flow equations have 
been used widely in the literature to study inflationary dynamics [11, 12, 13]. For 
a particular choice of initial conditions for the flow parameters, equations (4) can be 
integrated until the end of inflation without relying on any slow-roll assumption. In this 
way one can explore the full dynamics of a large number of simple inflationary models 
in a relatively model independent way. Nevertheless, in section 4.2, we shall explore 
another approach to inflation model building based on e(A^). 

Having generated a model of inflation, the spectrum of curvature perturbations 
can be calculated using analytic formulae such as those derived in [14, 15, 16, 17]. In 
this paper, we shall focus the accuracy of the Stcwart-Lyth formula [14], which has been 
widely used to calculate inflationary perturbations. This approach is not necessarily self- 
consistent if the conditions of slow roll are violated during inflation (which is, of course, 
allowed by equations (4)). An alternative approach is to calculate the perturbation 
spectra numerically using, for example, the Mukhanov formalism [18, 19]. We briefly 
compare these two approaches. 

2.2. The Stewart- Lyth approximation 

If the inflaton fleld (j) rolls so slowly that its kinetic energy is continually dominated 
by its potential energy, then the power spectrum of scalar curvature perturbations (we 
consider only scalar perturbations in this paper) can be approximated by the standard 
flrst-order slow-roll approximation: 

2 



m|i \H' 



(5) 

k=aH 
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where primes denote differentiation with respect to (f). The next order correction is given 
by Stewart and Lyth [14] 

C-3 \ 2 H"^ 



\H' 



(6) 

k=aH 

where C is the constant introduced in equation (3). Equations (5) and (6) are evaluated 
at the instant when each Fourier mode k crosses the Hubble radius. 

For numerical purpose, we express the Stewart-Lyth power spectrum purely in 
terms of the flow parameters: 

[l-e+^a^e-^ ( i'\n\k/ko\ \ 

Vn{k) ~ Vn{ko) ; ^ exp -2 / e(60 - N')dN' , (7) 

where ko = 0.002 Mpc^^ is the pivot scale at which the spectrum is normalised. This 
normalisation is estimated as 

Viziko = 0.002 Mpc"^) ^ 2.36 x 10"^ . (8) 

from the 3- year WMAP results combined with several other datasets [20]. 

2.3. The Mukhanov formalism 

The Stewart-Lyth formula assumes that the background parameters e, cr are 
approximately constant as each mode crosses the Hubble radius. However, as soon 
as the infiaton potential becomes steep, e and a will change rapidly and the formula 
becomes unreliable. To calculate the curvature spectrum accurately in such cases, we 
compute a numerical solution for the Mukhanov variable u [18, 19] 

a d(j) , . 

u = -zn, --H^, (9) 

where TZ is the curvature perturbation. In Fourier space, Uk evolves according to a 
Klein-Gordon equation with a time-varying effective mass: 

d'^Uk /, 9 1 d'^z\ , 

where r is the conformal time rfr = dtja. At early times r^, when the mode k is much 
smaller than the Hubble radius, Uk is usually chosen to represent the Bunch-Davies 
vacuum statej with 



During inflation, Uk evolves in time until the physical scale of each mode is far 
outside the Hubble radius, at which point \R\ approaches the asymptotic value and 
freezes. The power spectrum can then be evaluated as: 

k"^ Vi 2 



271^ 



k^aH 



X Some authors have suggested that trans-Planckian physics can be modelled by choosing a different 
vacuum state [21, 22] 
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To numerically integrate equation (10), we rewrite it in terms of the number of 
e-foldings N rather than conformal time 

2 

-fie,a,0 



- (1-e)— + 



dN^ 



dN 



k 

oH 



Uk = 0, 



(13) 



since we use the inflationary flow equations to calculate the evolution of the unperturbed 
background. The factor / in equation (13) is given by 



/(e,a,e) = 2 



, 3 

4e a 

2 



(14) 



and is exact despite its dependence on the flow variables. The initial conditions assigned 
to Uk are: 



1 



Re{uk{Ti)) 



Im{uk{Ti)) = , 



Re 



Im 



duk 
dN^ 

duk 

dN' 
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aH 



(15) 



where aH/k is a small fixed ratio, which we have set to be 10^'^. Wc have checked that 
the final power spectrum shows negligible change if this ratio is reduced further. 

For each inflation model, we continue evolving each mode until the end of inflation, 
so that 



Vn{k) = Vn{ko) 





3 

) 


Uk 


2 








end 



(16) 



As long as inflation satisfies slow- roll conditions (e, a -C 1), the numerical solution (16) 
will agree accurately with the Stewart-Lyth approximation (6) for all modes that are 
well outside the horizon at the end of inflation. However, there can be large differences 
if the slow conditions are violated temporarily during inflation. These types of models 
have been analysed by Leach et al. [23, 24]. In particular, it was shown in reference [24] 
that the Stewart-Lyth formula will fail if the quantity z in (9) becomes smaller than 
its value at the time a mode crosses the Hubble radius. Some specific examples of this 
behaviour are described in the next section. 



3. Limitations of slow-roll approximation 

As in previous papers [10, 11, 12] we use the inflationary flow equations (4) to evolve a 
large number of inflation models with initial flow parameters selected at random from 
the uniform distributions 

eo e [0, 0.8] , 
(70 e [-0.5, 0.5] , 

6 £[-0.05,0.05], (17) 
^^h\o e [-0.025 X 5-^+^ 0.025 x 5"^+^] , (3 < £ < 10) 
''A^lo = . 
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The distributions (17) should not be interpreted as a physical probability distribution for 
inflationary models, but merely as a device for generating models with a wide range of 
evolutionary histories. Despite the fact that the truncation is made at £ = 11, the types 
of models that we generated were sufficiently diverse for us to understand when and 
why the Stewart-Lyth approximation breaks down in general. Keeping the hierarchy to 
higher orders can produce more complicated inflationary models, but this will not affect 
our qualitative conclusions. 

We have selected three models for which the Stewart-Lyth approximation differs 
to varying degrees from integrating equation (10). First, we consider small scale 
perturbations where the Stewart-Lyth formula is inaccurate for modes for which \Tl\ 
does not have time to reach its asymptotic value well before inflation ends. Next, 
we describe an inflationary scenario where e changes rapidly, arising from a 'bumpy' 
potential. Finally, we show an example for which the potential has a steep feature, 
inducing fast roll and interrupting inflation temporarily. Models based on speciflc 
potentials V{(p) displaying similar behavior to each of these three models can be found 
in [23, 24, 4]. Throughout, we assume that there were 60 e-folds of inflation between 
the time CMB scale perturbations were generated and the end of inflation. 

3.1. Small-scale perturbations 

Figure la shows e{N) for a typical inflationary scenario that lasts for 60 e-folds generated 
from the distribution (17). Within the last few e-folds of inflation, e rises rapidly to unity. 
Consequently, the damping term in equation (13) is suddenly reduced, and Uk evolves 
more like a frictionlcss oscillator. This results in an upturn in the power spectrum Vn 
(flgure lb), peaking at the smallest scales that exit the Hubble radius. 

The numerical solution of the Mukhanov equation (10) (bold line) and the Stewart- 
Lyth formula arc almost indistinguishable for wavenumbcrs k < 10^^ Mpc^^ . Deviations 
in the solutions can be seen at larger wavenumbers, but they are relatively small. This 
behaviour is quite typical for models in which inflation ends when e exceeds unity. For 
most purposes the Stewart-Lyth formula will provide an excellent approximation to the 
power spectrum even on small scales. Large deviations from the Stewart-Lyth formula 
require more unusual inflationary models such as the two examples described next. 

3.2. 'Bumpy' potential 

Figure 2a illustrates a scenario in which e peaks strongly during inflation as result of a 
potential with a pronounced bump or drop. Here the Stewart-Lyth approximation and 
the exact integration agree at large scales but differ considerably as the background 
evolves rapidly leading to the prominent dip in V-r, as seen in figure 2b. We can 
understand this feature as follows. As the background value of e rises around a 
potential bump during infiation, the dynamics of the fiow variables require that a 
simultaneously plummets to large negative values. Thus, according to the slow-roll 
formula for the spectral index (3), the gradient of the power spectrum also decreases 
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Figure 1. Panel (a) shows the parameter e as a function of N, the number of e- folds 
of inflation (with A'' = at the end of inflation). Rise in e towards unity at the end of 
inflation gives rise to the small scale feature of the power spectrum Vnik), shown in 
panel (b). The spectrum in bold line uses the Mukhanov variable (16) while the thin 
line shows the Stewart-Lyth approximation (7). 




Figure 2. Panel (a) shows e{N) for a model in which e(iV) peaks prominently during 
inflation as a consequence of a bumpy potential with a pronounced bump or drop. The 
origin of the dip in the Mukhanov spectrum (panel (b), bold line) and its magnitude 
in comparison with that of the Stewart-Lyth spectrum are discussed in the text. 



to large negative values. This explains the downturn in the Stewart-Lyth spectrum. 
However, superhorizon evolution, which is neglected by the Stewart-Lyth approximation, 
becomes important when the background changes rapidly. In this a becomes 

large and negative, one finds that the Mukhanov variable evolves as ~ e"'"^'^ where 
A = 0{a). This results in an additional decrease in Vn around the bump. Subsequently, 
towards the true end of inflation, the tail of the spectrum resembles that in flgure lb as 
expected. 
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Figure 3. Panel (a) shows e(N) for a model which temporarily fast-rolls (when e > 1), 
giving rise to a typical 'ringing' in the power spectrum shown in panel (b). The origin 
of the large discrepancy between the Mukhanov (bold) and the Stewart-Lyth spectra 
in panel (b) is discussed in the text. 



3.3. Fast roll 

Consider a scenario in which e(A'^) peaks so strongly during inflation that e exceeds unity 
temporarily, hence breaking up inflation into two or more stages (see e.g. [25]). Such a 
'fast-roll' model is shown in flgure 3. 

The power spectrum resulting from fast roll (figure 3b) shows a large difference 
between the Mukhanov and the Stewart-Lyth spectra. In addition to the enhanced dip 
as explained earlier, wc draw attention to two other features: i) the subsequent sharp 
rise of the Mukhanov spectrum several orders of magnitude above the Stewart-Lyth 
approximation, and ii) the rapid oscillating features across 10^^ < /c/Mpc^^ < 10^^. 

These features can be understood as follows. Equation (13) clearly shows what 
happens to the Mukhanov variable Uk during fast roll. The (1 — e) friction term 
of the harmonic oscillator becomes negative during fast roll, hence boosting \uk\ 
momentarily. Physically, fast-rolling leads to an entropy perturbation which sources 
growth of curvature perturbations on superhorizon scales [23]. The finer-scale 'ringing' 
induced during fast roll have previously been seen in [6, 4] in their investigations of a 
potential with a step. 

4. Application to primordial black holes 

Many authors have speculated on the possible existence of primordial black holes 
generated during the early Universe [26, 27, 28, 29]. Such objects, if they are more 
massive than 10^^ g could contribute to the dark matter density at the present time 
[30, 31] and may be the source of high energy cosmic rays [32, 33] (see [34] for a recent 
review). However, it has proved difficult to think of compelling scenarios in which an 
interesting density of primordial black holes could form. Formation of primordial black 
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holes in inflationary models have been discussed by a number of authors (sec [35, 36] and 
especially [37], which is the closest in spirit to this work). We revisit this problem in this 
section and show that significant primordial black hole formation during infiation could 
have occured only if the infiationary dynamics were highly unusual and finely tuned. For 
such unusual inflationary models, an accurate computation of the primordial black hole 
abundance requires numerical integration of the inflationary perturbations as described 
in the previous section. 

4.1. Basic quantities 

Here, we consider black holes that form when inflation ends and primordial density 
fluctuations begin to re-enter the Hubble radius. In this section we review the basic 
quantities relevant to the formation of primordial black holes and compute their relic 
abundance using the Sheth-Tormen mass function [38]. 

Fluctuations in the energy density smoothed over a spherical region of scale R are 
related to the energy density fleld S{x) by the convolution 

S{R,x) = y"l^(|x'-x|/i?)(5(x')ciV, (18) 

where W denotes a 'window function' of scale R . In Fourier space, the variance of 
6{R, x) is given by 

Vsik)W'{kR)^, (19) 

where W{kR) is the Fourier transform of the window function W. Here we will assume 
that is a top-hat function of radius R (see e.g. [39], §4.3.3). 

The power spectrum of the energy density perturbations Vs{k) = A;^/27r^(|5(/c)p) 
is related to the spectrum of curvature perturbations Vn{k) by [39]: 




A region with density contrast 6 forms a black hole if S exceeds a critical value Sc- 
Previous authors have shown that 5c should be roughly equal to the equation of state 
during radiation era {Sc — 1/3) [28, 40, 41], and we shall take Sc — 1/3 throughout this 
paper. § In addition to the lower limit there is an upper bound for 5 of approximately 
unity, since for such high overdensities large positive curvatures will cause patches to 
close on themselves, creating disconnected 'baby' universes [28]. 

Next, to obtain the mass spectrum of primordial black holes, consider the mass M 
associated with a region of comoving size R and average energy density p. Early work 

§ but see [42, 43, 44] for differerent opinions on the value of 6c inferred from simulations of gravitational 
collapse. 
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by Carr [45] suggested that primordial black holes form with mass comparable to the 
'horizon mass', i.e. 

M^—w'/'p{H-')\ (21) 
3 

Recent numerical investigations [42], however , suggest that M may be a small fraction 
of the horizon mass, and may be bounded from below [44, 46]. Nevertheless, these 
results have been criticised by some authors [43], and we shall still use equation (21) in 
this paper. It then follows that M and R are related by [47] 



R / M\^^'^ /n \ 1/12 

^ -2.2x10-2^ — 1 



1 Mpc V 1 g / 

where the effective relativistic degree of freedom is of order 100 in the very early 
universe when the temperature exceeds ~ 300 GeV. By combining Equations (19)-(22), 
the variance cr{R) can then be converted to the mass variance a{M). 

The fraction of the universe, ilpbh, in primordial black holes of mass greater than 
M is 

Ppbh 



pbh 



[ P{5{M))d5{M), (23) 

Ptotal JSc 

where P{S) is the density fluctuation probabihty distribution smoothed on scale M and 
is such a sharply falling function that the upper limit can be extended to infinity without 
loss of accuracy. Previous work [48, 49, 50] has used the Press-Schechter formalism, in 
which P{5) is a Gaussian function, yielding 

Here we use the Sheth-Tormen mass function [38] which is based on an ellipsoidal model 
of gravitational collapse and provides a better fit than the Press-Schechter function to 
the abundances of non- linear structures in numerical simulations [51]. The primordial 
black hole abundance now becomes: 



f^pbh(> M)=A 



\ , 1^/1 aS', \ 



(25) 



where A ~ 0.3222, p ~ 0.3, a — 0.707 and r(x, y) is the incomplete gamma function. 



4.2. Results 

For inflationary models that support primordial black hole production, we calculated 
the black hole abundance using the Sheth-Tormen formula] | (25) with M — lO^^g, 

II In calculating ilpbhj the Sheth-Tormen and the Press-Schecter formalisms broadly agree. However, 
this depends sensitively on the shape of the inflaton potential. We found classes of potentials where 
the two prescriptions differ by a few orders of magnitude. Nevertheless, the resulting bound on ris (30) 
remains valid regardless of which prescription is used. 
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corresponding to the lightest black holes that may survive till present day without 
evaporating completely. For models that are radiation dominated at the end of inflation 
(i.e. do not experience an extended matter dominated phase), the observational 
constraints on the primordial black hole abundance can be described roughly by 

npbh(M > W g) < 10-20. (26) 

More accurately, the constraint on flphh is a function of the mass of the primordial black 
holes and on the physics of the (p)rchcating phase [52, 53]. In this paper it will be 
adequate to use the simple criterion (26) to give a feel for the types of inflation models 
that might produce black holes. 

4.2.1. Models with blue spectra A viable model of inflation requires at least 60 e- folds 
of inflation [54]. Using the inflationary flow equations, we flnd that the vast majority 
of inflationary flow models that experience such a large number of e-folds are models 
in which the inflaton is trapped in a local minimum of the potential [10, 11]. For such 
models, inflation never ends unless some additional physical mechanism is included. 
Examples of such mechanisms include an auxiliary field, as in hybrid infiation [55] or 
open strings becoming tachyonic in brane infiation models [56]. In these models, the 
resulting spectrum of curvature perturbation is very well approximated by a power- law 
relation: 

Vnik) = Vniko) f , (27) 



\ko/ 

with spectral index ris given by equation (3). As the inflaton sits in the local minimum, 
e remains very small while a is positive and may be large. Thus, the power spectrum 
for this type of models is 'blue' (ris > 1). The running of ris is given by 

dns 1 2 r /^3 - 5C\ ^ /3 - C\ \ 

'2i- 28e^ -5ea-l + —— ) ] . (28) 



dink 1-e \ \ 2 J ^ \ A ^ ^ 

Because the leading order comes from the small ^ term, the running in these blue models 
is negligible. Therefore, the shape of the curvature power spectrum in these models can 
be characterized accurately by a constant spectral index Ug. 

Consequently, for models with blue spectra, Ug directly determines the abundance 
of primordial black holes Qpbh, as shown in Figure 4. The points show models trapped 
in a local minimum of the potential for which inflation was arbitrarily terminated after 
200 e-folds. The curve shows the fltting function 

logio ^Pbh = -«(! + {(3{ns - 1.295))0"'^'' , (29) 

where a ^ 30, (3 = 10.2891, p = 1.05307 and q = 0.218219. Thus, the constraint 
^pbh ^ 10-2° discussed above implies an upper bound in ris 

ris < 1.3, (30) 

in excellent agreement with the analysis of [52] , which used constraints on the amplitude 
of the scalar fluctuations derived from the four-year COBE data [57]. However, the 
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Figure 4. Plot of the primordial black hole abundance Jlpi5h(Af > lO^^g) against the 
spectral index Us for 10^ models where inflation ends suddenly. There is virtually no 
dispersion around the curve given approximately by equation (29). The exponential 
sensitivity of ilphh "-s at the steep part of the curve highlights the level of fine-tuning 
needed to produce an interesting density of primordial black hole . 

observational constraints from the 3-year WMAP and SDSS [20, 58, 59, 60] can now 
convincingly exclude a high density of primordial black holes from this type of model 
(for example, reference [58] quotes = 0.964 with a Icr error of only ±0.012). 

4-2.2. Designing potentials with t(N) The abundance of primordial black holes is 
closely linked to the power spectrum Vn{k) on small scales. If the power spectrum 
has a high amplitude at small scales, as in the case of blue spectra described above, 
a significant density of massive primordial black holes may form. However, using the 
inflationary flow approach, models that support 60 e-folds of inflation and end with 
e = 1, are in general, tilted 'red' over a large range of scales and therefore cannot 
produce primordial black holes. Even the small upturns seen in figures lb and 2b occur 
at such low amplitudes that fipbh remains negligible. 

However, as discussed in §3.3, for models which experience a period of fast roll, it is 
possible to generate a large jump in Vn at small scales, and this can lead to primordial 
black hole production. For such models, it is important to use a numerical computation 
to calculate the power spectrum rather than the Stewart-Lyth approximation. For 
example, for the model shown in figure 3, the Stewart-Lyth approximation predicts 
a negligible density of primordial black holes whereas integration of the Mukhanov 
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Figure 5. Left panel shows e(A'^) designed using Hermite interpolating polynomials. 
Fast-roll regime has been modelled so as to produce primordial black holes at an 
interesting abundance. The reconstructed power spectrum Vn{k) is shown in the 

right panel. For this model, fipbij(Af > 10^^ g) ~ 1.6 x 10"^". The observables at 
CMB scales are r ~ 0.01, ~ 0.99 and dug/dhik ~ 10~^, broadly compatible with 
constraints from WMAP and other data. 



equation gives Qpbh(-/W^ > lO^^g') ~ 10~^. One can therefore ask if it is possible 
to construct fast-roll models that can produce an interesting primordial black hole 
abundance, while satisfying the current observational constraints on the shape of the 
power spectrum at large scales. 

We can construct such models easily based on the Hamilton- Jacobi approach by 
specifying the dynamics of inflation in terms of e{N) ^. We first choose a set of values 
of e{N) with the boundary conditions: 



e(0) - 1, lim e = 0+. (31) 

AT— Kx) 

The observational constraints on r, ris and its running can be satisfied by adjusting 
the values of e, de/dN and d'^e/dN'^ a,t N — 60. The form of e(A'") can then be adjusted 
'by hand' to generate a brief period of fast roll at small values of N to produce a 
significant density of primordial black holes. 

We then construct e(A^) by using piecewise cubic Hermite interpolating polynomials. 
These polynomials are well suited for our purposes since it is easy to apply constraints 
on the positions of cxtrcma and on the derivatives. The values of e{N) and its first and 
second derivatives then completely determine the evolution of the Mukhanov variable 
u via equation (13), and thus the power spectrum Vn can be reconstructed. Figure 5 
shows e(A'") for such a fast-roll model and the reconstructed power spectrum. For this 
model, 

Qpbh ^ 1.6 X 10-2°, r~0.01, 71^-0.99, dris/dlnk lO'^, (32) 

% Bond J.R., private communication 
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and so is broadly compatible with observational constraints [20] yet produces a 
significant abundance of primordial black hole . It is, of course, possible to generate 
many variants of this type of model with different values of and dus/dlnk on CMB 
scales. But all such models will require a finely tuned period of fast roll towards the 
end of infiation if they are to produce remnant black holes. 

We end this section with a brief comment on the effect of non-gaussianities on 
primordial black hole abundance. Chen et al. [61] recently pointed out that step-like 
features in the inflaton potential can give rise to large non-gaussianities in the primordial 
perturbations. The effect of non-gaussianities on the abundance of primordial black 
holes was considered in [62, 63]. These authors found that+r 



where Q — C(/nl)- Thus, even large non-gaussianities anticipated in fast-roll models 
will be exponentially suppressed, and the overall effect on flphh is expected to be small. 

5. Conclusions 

There have been many papers discussing the formation of primordial black holes from 
inflation. The simplest such models involve slow-roll inflation and inflationary dynamics 
that lead to a blue perturbation spectrum. For this type of model, our results conflrm 
earlier work that shows that the scalar spectral index must be less than rig ~ 1.3 to 
avoid overproduction of primordial black holes. However, we also argue that the most 
natural way to realise such models is to arrange the infiaton to be trapped in a minimum 
of a potential and for infiation to end abruptly, as in hybrid infiationary models. For 
such model, the run in the spectral index and the amplitude of the tensor component is 
negligible. Therefore, observational constraints on the scalar spectral index determined 
on CMB scales can be extrapolated to the small scales relevant for primordial black 
hole production. The strong constraints on the scalar perturbations from WMAP and 
other data therefore rule out primordial black hole production from slow- roll inflationary 
models. 

It is, however, possible to generate infiationary models which satisfy observational 
constraints on CMB scales yet produce a significant density of primordial black holes 
provided there is a period where slow-roll conditions are violated towards the end of 
infiation. We have shown in section 3 that for this type of model, it is essential 
to compute the perturbation power spectra numerically, rather than using slow-roll 
approximations. Section 4.2 presents a construction method for such a model based on 
designing e(A^). While it is possible to produce primordial black holes from inflation, 
the dynamics of this type of fast-roll model must be very flnely-tuned to produce an 
acceptable density of primordial black holes. In general, even if a model is constructed 
to have a period of fast roll, the primordial black hole density will either be excessively 

+ J.-C. Hidalgo, private communication. 




(33) 
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high or neghgible. Although there has been renewed interest in primordial black holes 
as a possible explanation for the existence of intermediate mass black holes [64], and 
to account for the rapid growth of supermassive black holes in the centres of galaxies 
[65] , it seems extremely unlikely to us that primordial black holes formed as a result of 
inflationary dynamics. 
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